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Q_ij Abstract 

We investigate charm production in an equilibrated quark-gluon plasma pro- 
duced in heavy-ion collisions at RHIC and LHC. Effective quark and gluon masses 
are introduced from thermal QCD calculations. Assuming a Bj or ken- type longitu- 
^ i dinal expansion and including the influence of temperature dependent masses on 

the expansion, we determine the total number of cc pairs produced in the quark- 
gluon plasma phase. We calculate the charm production rate at leading order with 
massive gluons and quarks and compare our result to charm production by massless 
partons. We consider two different scenarios for the initial conditions, a parton gas 
with a rather long kinetic equilibration time and a minijet gas with a short equi- 
libration time. In a parton gas, assuming m c = 1.2 GeV, we obtain a substantial 
enhancement over the thermal charm rate from massless quarks and gluons, up to 
4.9 secondary charm quark pairs in Au+Au collisions at RHIC and 245 charm pairs 
in Pb+Pb collisions at the LHC. 
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1 INTRODUCTION 

Two decades ago, studies of quantum chromodynamics (QCD) predicted that quarks 
and gluons can be liberated from hadrons and, at appropriately high energy densities, 
a new state of the matter, the quark-gluon plasma (QGP) will be formed Lattice- 
QCD calculations have shown that the critical energy density of this phase transition is 
e c ~ 2 — 3 GeV/fm 3 . The critical temperature is T c m 260 MeV for a pure gluon plasma 
while T c « 150 MeV when dynamical light quarks are included 0. These conditions are 
expected to be produced in ultrarelativistic nucleus-nucleus collisions. Recent experiments 
have not unambiguously proved QGP formation at the CERN SPS where = 17-20 A 
GeV, but in heavy-ion collisions at RHIC (y/s = 200 A GeV) and LHC (y/s =5.5 A TeV) 
the required energy densities should be produced. 

At ultrarelativistic energies, nucleus-nucleus collisions can be described by parton 
interactions in perturbative QCD inspired models ||, f|]. In this framework, hard and 
semi-hard scatterings among partons dominate the reaction dynamics. These collisions 
are able to drive the parton system toward thermal equilibrium very effectively || EJ, M], 
but chemical equilibrium may not be established || |7|]. Equilibration strongly depends 
on the actual parton densities in the colliding nuclei and the microscopic processes on 
the parton level. The behavior of the parton densities at small x is not fully established: 
newer parametrizations of the parton densities already greatly increase the energy density 
of the initial state M and the introduction of softer semi-hard scatterings (a smaller p? 
cutoff on the minijet cross section) could yield a parton system very close to equilibrium 
a short time after the primary collisions. 

During the time-evolution of the parton system, in parallel with the thermalization, 
many collective phenomena can appear. The massless color degrees of freedom, the quarks 
and gluons, are strongly interacting, as described in the framework of finite temperature 
QCD ||. In this theory, temperature dependent thermal quark and gluon masses are 
introduced which can dramatically affect the phenomenological description of the time 
evolution. Massive fermionic and bosonic excitations can be derived in many ways, e.g. 
by considering leading order medium effects on the QCD self-energy || . Since the behavior 
of the QGP and its high-temperature excitations are not well understood, it is challenging 
to find measurable tests of this behavior in heavy ion collisions. 



In general, open charm [JT0], [IT], [12], [13], |14| , direct photons and dilepton production 
15| can be used as direct probes of the early parton dynamics and the evolution of 
the QGP. Primary charm production from the initial nucleon-nucleon collisions has been 
calculated and can be substantial |]TTJ , especially when up-to-date parameterizations of the 



parton densities which increase at low x, are used [11, 13] . In the most recent calculation 



of initial charm production |13|] in Au+Au collisions, extrapolated from pp collisions, a 



total of 9 (450) initial cc pairs were produced at RHIC (LHC), primarily due to the 
behavior of the parton densities at low x. As shown in ref . [|13| , shadowing can reduce the 
yield in the central region by up to a factor of two, depending on the cc pair mass and 
transverse momentum. 
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The amount of secondary charm production in the QGP phase is an open question 
fT2| , p~3| , |14| , p!6| . If the parton system is dilute || [L4j], secondary charm production would 
be suppressed because of the low parton fugacities. Furthermore, at RHIC the charm mass 
is 2-3 times larger than the expected temperature scale, T ~ 0.35 — 0.55 GeV, and much 
greater than the bare light quark masses. Thus in this case, even in a fully equilibrated 
plasma charm production is not significantly enhanced. At the higher LHC energy, the 
predicted initial temperatures, T « 0.7 — 1 GeV, are close to the charm mass. One 
can then expect a larger thermal charm production rate than at lower energies, perhaps 
comparable to the initial yield. 

Charm production in the plasma by massless quarks and gluons is small because 
the massless partons must be very energetic to overcome the cc pair mass threshold. 
However, the effective thermal quark and gluon masses generated by the plasma could 
significantly enhance the thermal charm yield, perhaps even approaching the initial yield 
after shadowing has been included. This additional charm production is a possible probe 
of massive excitations in the QGP which could test the level of thermalization and the 
evolution of the quark gluon plasma. Furthermore a large charm quark multiplicity would 
favor the production of exotic phenomena such as multi-charm hadron production p 



The influence of massive gluons on strangeness production has already been inves- 
tigated [Tj| and found to have a negligible influence on the total strangeness production 



rate. However, that calculation considered massless light quarks and a temperature re- 
gion, T w 200 — 250 MeV, only slightly above T c . Here we perform a similar calculation 
but with massive quarks as well as massive gluons and with the higher initial tempera- 
tures obtained from more recent estimates of the initial conditions J7|, Under these 



circumstances, the thermal quark and gluon masses can generate a significant increase in 
secondary charm production. 

In this paper we investigate the thermal charm production rate in a fully equilibrated 
QGP with massive quark and gluon degrees of freedom. We compare our results with 
those from zero mass quarks and gluons. We will assume both a pure gluon gas and 
a quark-gluon plasma. We also distinguish between parton gas and minijet gas initial 
conditions. We obtain the total number of cc pairs produced during the lifetime of the 
plasma by assuming a Bjorken-type longitudinal expansion for Au+Au collisions at RHIC 
and Pb+Pb collisions at LHC. The influence of temperature dependent masses on the 
expansion and the speed of sound is also considered in detail. 



2 EFFECTIVE QUARK AND GLUON MASSES 

In a strongly-interacting QGP, collective fermionic and bosonic excitations appear, as 
also seen in QED. These quasi-particles can be interpreted as massive quarks and massive 
gluons. The propagators are modified in the medium because the color charges are dressed 
by the interaction with their neighbors in the plasma phase. The effective masses are 
related to the properties of the QGP and can thus be characterized by their temperature 
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and density dependence. Using the high temperature approximation [2C] or, equivalently. 



the hard thermal loop approximation |2T| which ensures a gauge invariant treatment, one 
can obtain the following effective equilibrium masses at zero chemical potential, /i = 0, 
for quarks and gluons massless at T = |§: 



o 2 T 2 

n 2 T 2 AT 

< th (T) = + (2) 

We assume that the equilibrated QGP contains light u, d, and s quarks and anti-quarks 
so that Nf = 3 and m gt th = gT / Since the strange quark has a non-zero bare mass, 
we consider an approximate effective mass, rh s , where 



™> s = V m ^o + m g,th > ( 3 ) 

and m Sj o = 150 MeV. We neglect the small current mass of the light quarks so that 
rh q m 9j th- The effective mass of the gluon is similar, rh g = m g ^h- Charm quarks are 
assumed to be produced with their bare mass and dressed later. The hadronization of 
charm quarks does not influence the charm production rate. 

The temperature-dependent running coupling constant is 

9 KtJ (33 -2iV / )ln[(19T c /A M s) (T/T c )\ ' 1 ] 

where in SU(3) gauge theory, T c /Ay^ = 1.78 ± 0.03 f22fl . With this dependence, the 
thermal gluon and quark masses increase approximately linearly with in the temperature, 
m g ,tb ~ (130 T/T c + 36) MeV and m qjth w m g ^/V^, where T c = 150 MeV. Then at 
T w 500 MeV, ih g w 470 MeV while at the lower temperature assumed earlier fl8fl , 



T ~ 300 MeV, rh g ~ 300 MeV. In ref. |TS|, 53] massive gluons were introduced by 



reproducing lattice-QCD energy density and pressure results [p^j . The effective gluon 
mass, assumed constant, was M g ps 450 - 500 MeV for 1.2 < T/T c < 2.4, 30 - 40% larger 
than our value in the same region. 



3 CHARM QUARK PRODUCTION 

The processes relevant for cc pair production in the quark-gluon plasma are the following: 

q + q — > c + c (5) 
s + s — >■ c + c (6) 

g + g — > c + c (7) 

An analogy can be drawn between charm and strangeness production by the plasma. In 
early calculations of strangeness production, processes (5) and (7) were treated at leading 
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order with bare quark and gluon masses [Q. On the other hand, effective quark masses 
were used as infrared cut-offs simultaneously removing electric and magnetic infrared 
singularities, see e.g. ref. ||. Effective gluon masses were included in calculations of 
chemical equilibration [g] and strangeness production fl8j in the QGP, and in calculations 
of radiative energy loss in a parton gas [26| . In these cases the bare Green's functions were 
used with the effective gluon mass but the quarks were not dressed. Here we combine 
these approximations, using the bare Green's functions containing dressed quarks and 
gluons and calculating the production rates in perturbative QCD at leading order. 



In the calculation of strangeness production with massive gluons one additional 
process was included, the decay of massive gluons. From the thermal gluon mass, eq. (0), 
and the temperature dependent running coupling constant, eq. (f|), we estimate that the 
threshold for g — > cc, rh g > 2m c , will be reached only at T « 20 T c . We therefore neglect 
this channel. We can also assume that the annihilation of produced cc pairs is negligible. 
The total charm production rate is thus: 



dN cl 
dr 



[A g (T) + 2A q (T) + A s (T)].V(r) 



where 



MT) = — 
nA 



ds / 2aAs)\ sis — Am 2 



ME,) fj(E 2 ) d 3 Pl d 3 p 2 



(8) 



(9) 



2Ex 2E 2 (27r) 3 (27r) 3 

and Hi is the number of identical particles in the initial state. We calculate A q and A s 
separately since m q ^ m s . The Bose distribution is used for gluons, f g (E) = [e E l T — l) -1 , 
and the Fermi distribution is used for quarks, f g ,q(E) = (e E ' T + l) -1 . We consider a 
symmetric QGP with zero baryon and strangeness chemical potential, fi q = fi s = 0. The 
total production cross section for each channel, a,, is calculated at center-of-mass energy 
y/5. 

The gluon-gluon fusion and quark-antiquark annihilation rates can be written in the 
Lorentz-invariant form 

MT) = ^Jdsf^f*ZL2 a(s) Js(s - 4m?) S(s - ( Pl +p 2 ) 2 ) 



(2tt) 3 7 (2tt) 3 
•8(pi-rnZ)5(jZ-mZ)MEi)ME2) 



(10) 



with four momenta p\ = (Ei,pl) and p 2 = (E 2 ,p%). Following ref. |T8[ this integral can 
be rewritten as 



A(T) = -— — / dsa^sis-Am^Y: E(±! 

6Z7T* TlA JAmJ ,~ r.r,n 



oo 1-2 



2 k=-l+2 



'IF 



where the sum over index k is incremented by 2, with the + sign used for bosons and the 
— sign for fermions. The modified imaginary Bessel function, Ki(a,ki), has the argument 



Q-kl 



2T\ 



k 2 



(12) 
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The use of a simple Boltzmann distribution instead of the Bose and Fermi distributions 
is equivalent to reducing the sum to only the I = 2, k = term with a 02 = y/s/T. At high 
temperatures and high effective quark and gluon masses, the Boltzmann distribution is a 
good approximation. However, here we will use the full sum. 

Once the time evolution of the volume, V(r), and the temperature, T(r), are de- 
termined, we can calculate the appropriate quark and gluon masses and the total charm 
pair production rate. 



3.1 ELEMENTARY CROSS SECTIONS 

Here we explain the calculation of the total cc subprocess production cross sections for 
massive quarks and gluons. As discussed earlier, we neglect cc pair annihilation as well 
as the decay of massive gluons. 

When the light quarks have an effective mass, the square of the matrix element for 
processes (4) and (5) is 

IM-I2-9 4 Tr [7 M (^-^)7 t '(^ + ^)] Tr[7 M (j/ 1 -m c )7 ;y (j/ 2 +m c )] 

1 ll ~ 9 (s-mi + ryiy + rimi ' 1 } 

where p^ and are the four-momenta of the incoming quarks, % = u, d and s, and pi 
and P2 are the four-momenta of the outgoing c and c quarks. The propagator has been 



modified by the finite mass and width of the gluons as in ref. fL8fl . After evaluation of 
the traces, the square of the matrix element can be expressed as 

\M\* ,,2^2 64 (m\ + mf - t) 2 + (m 2 c + mj - u) 2 + (2m c 2 + 2mj)s 

Ivh = a ; ir a fl ; — —r-. . _ — ^ n - , 14 

11 3 9 (s -m2 + r 2 /4) 2 + Y 2 m 2 g ' v ; 

where g 2 = 47ras and d% = 6 for the spin and color degrees of freedom for each flavor. 
Note that eq. ([13]) does not include the usual average over the initial spin and color while 



eq. (|T4l) does include these factors. 

The total cross section is obtained from the t integration of |Mj| 2 



CTi(S) 



1Qtts(s — 4mj) 
with the limits 



dt \Mi\ 2 (15) 




where the sum of the Mandelstam invariants s, t and u is 

s + t + u = 2m 2 + 2m 2 c . (17) 
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For gluon fusion, (7), the total cross section can be obtained from the integral 

1 /■* 



g '<' )= lfrr.(,-4rf>j) /- + + < 18 > 



with the integration limits 



/s — 4m 2 , / s _4 m 2\ 
Detailed calculations of the invariant matrix elements, M s , M u and M t , are given in ref. 



T8|j for ss production by transverse gluons of constant mass. The same expressions can 



be used here with the change m s — > m c . The sum of the Mandelstam invariants is now 

s + t + u = 2m 2 , + 2m 2 . (20) 



We use two different values for the charm quark mass: m c — 1.2 GeV and m c 
1.5 GeV. The lower value was found to produce agreement between charm production 



calculated to NLO O] and pp total cross section data. This value has also been used in 



recent estimates of the charm contribution to the dilepton yield [I3[ . The larger value is 



somewhat more standard [|T|, [10, [L4| and, at energies near the charm production threshold, 
allows an all-order resummation of soft and virtual gluon corrections p7|. For smaller 
charm quark masses, the series cannot be resummed. 



3.2 THE HYDRODYNAMICAL MODEL 

We now discuss our calculations of the equation of state and the time evolution of the 
plasma. The pressure of an ideal gas of massive particles is 

gi r°° dk k 4 



. 6tt 2 Jo Jk 2 + m 2 



p = Y,-r-J i n [exp^v/F + m^Ti]- 1 , (21; 



where i — g,u, u, d, d, s, s. In the case of a symmetric plasma, /i n = fid = = 0, and we 
have Pi = Pj. Since rhi is temperature dependent, the energy density, e = TdP/dT — P, 
has an additional term proportional to the mass gradient 



9i I °° dk k r t a | ~Ij2\ _ ii— 1 ( i,2 i _i-_2 rn^^i 

. 27r 2 io JW+mt 

We assume a simple longitudinal expansion and follow the time evolution through 
the entropy 

s(t) = s(to)(t /t) • (23) 
We also calculate the square of the sound speed, defined as 

c 2 = ^ (24) 
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where dP/dT = s = (e + P)/T as a check on how closely the evolution follows that of an 
ideal massless gas with c 2 = 1/3. The temperature gradient of the energy density is 



de 
dT 



E 



1 9i 
T2vr 2 



dk 



k 2 + m? 



[exp ((3Jk 2 + m 2 )Tir 1 {(k 2 + 



m 



2\2 



+ (k 2 + m 2 ) 3k 2 - 2m { T 



drhi 
~dT 



(25) 



T 



[rut—) I [~dT 



A d 2 rhi rhi drhi 
+ m i —— r + 2- 



c/T 2 



T dT 



Note that the above expressions reduce to c 2 = 1/3 when rhi is constant. 

The mass gradients can be expressed rather simply when rhi — cg(T)T and the 
current quarks and gluons are massless. In this case we have 



drhi 

It 

d 2 rhi 
HT 2 



mi 
T 



1 



1 (rhi 

2 [dT 



1 rh- 



2c 2 T 4 



1 



3 (rhj 
2 \cT 



(26) 
(27) 



For the strange quark, rh s = yml + m 2 (T), the derivatives are 

dm* 



1 dm n 



dT 

d 2 rh s 
dT 2 



rh* dT 



1 / dm,, 



m 3 „ \ dT 



+ 



1 d 2 m q 
m, dT 2 



(28) 
(29) 



where dm q /dT = drhi/dT as above. Including effective quark and gluon masses tends to 
slow the evolution of the system as well as increase the cc rate. 

We need to fix the space-time volume to calculate the absolute number of produced 
charm pairs. The particle number is obtained from 



N = J ri^do^ = pirR 2 r J drj cosh rj = pV 



(30) 



where V = 2ir R 2 t sinh^*. To determine the maximum space-time extent of the plasma 
in rapidity as a function of r, we use the total available energy as a rough estimate: 



tot 



2AE] 



beam 



T "da,, 



(31) 



/rf 
dr][e cosh 2 r] + P sinh 2 rf\ . 
-n* 



From a comparison of the calculated E tot with the available energy, 2AE^ eam , we can 
determine the value of rj* and calculate the volume. With this volume, the charm yield 
is calculated from eq. (H). 
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4 TIME EVOLUTION OF THE SYSTEM 



We will consider the hydrodynamical evolution of a fully equilibrated plasma with two 
different sets of initial conditions, T and r , for RHIC and LHC collisions. 

The first parameter set is based on the parton gas model derived from the HIJING 
Monte Carlo code |§. The kinetic equilibration time, relatively long, is reached when the 
momentum distributions are locally isotropic due to elastic scatterings and the expansion 
of the system, r ~ 0.5 — 0.7 fm. We use T = 550 MeV and r = 0.7 fm at RHIC and 
T = 820 MeV, r = 0.5 fm at LHC, as in the ideal case described in ref. []T9j . With 
this model we do not distinguish between a gluon gas and a quark-gluon plasma in the 
estimate of T . 

The second set of initial conditions was determined from estimates of minijet pro- 
duction 0. In this case, the momentum scale of the minijet calculation sets the initial 
time, ro ~ 1/pt < l/po ; and the minijet yield determines the initial temperature. For 
a typical value of the momentum scale, po ~ 2 GeV, ro ~ 0.1 fm. Because the minijet 
yield depends on the composition of the system, the initial temperature depends on the 
partonic degrees of freedom. Since minijets predominantly produce gluons, if we consider 
only a gluon gas, the highest T is obtained, T = 445 MeV at RHIC and 1140 MeV 
at LHC. When light quark production is included, the number density increases but T 
decreases to 360 MeV at RHIC and 900 MeV at LHC. The plasma resulting from the 
early equilibration time and high temperature is referred to hereafter as a minijet gas. 

We compare and contrast the evolution of the plasma with the parton gas and the 
minijet gas in Au+Au collisions at RHIC in fig. 1 and in Pb+Pb collisions at the LHC in 
fig. 2. In each case, we show results for a pure gluon system and a quark-gluon system 
both with massless and massive partons. The time evolution of the temperature, energy 
density, square of the sound speed in the medium and the volume of the plasma are given 
in each case. 

The parton gas has a longer lifetime due to the longer equilibration time. The 
temperature evolution is shown in figs. 1(a) and 2(a) for the parton gas and in 1(b) and 
2(b) for the minijets. The difference between the temperature evolution with gluons alone 
and for the quark-gluon gas is seen to be small, independent of whether or not the partons 
are massive. The evolution slows when the parton masses are finite. We have cut off the 
evolution at T c = 150 MeV. The parton gas, with its slower evolution, remains above 
T c for t < 10 fm. At RHIC energies the temperature of the minijet gas drops below T c 
at r ~ 1.5 fm in the quark-gluon system. The pure gluon system, with its larger To, 
remains above T c for at least 2.5 fm at RHIC. For both the parton and minijet gas, the 
higher initial temperatures at the LHC make the finite parton masses more effective in 
slowing the evolution, in part because the effective parton masses are larger for the higher 
temperature. 

The increase in the number of degrees of freedom between a pure gluon gas and a 
quark-gluon gas is clearly reflected in the difference between the energy densities of the 
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two systems, shown in figs, l(c-d) and 2(c-d). Note that in both cases at RHIC, the 
energy density of the massive gluon gas is actually reduced relative to the massless gluon 
gas. This is due to the temperature gradient of the gluon mass. A similar effect is also 
observable in the minijet quark-gluon gas. In all these cases, the temperature is not large 
compared to the effective masses, causing the reduction. At the much higher temperatures 
of the LHC, although the effective masses are increased, the energy density is increased by 
the inclusion of the finite parton masses. We note that at RHIC, the energy density of the 
minijet gas actually drops below 1 GeV/fm 3 for r > 1 fm, even though the temperature 
remains above T c until r ~ 1.5 fm in the quark-gluon system and until r ~ 2.5 fm in the 
gluon gas. This suggests that the assumption of an equilibrated minijet gas at RHIC is 
perhaps questionable. 

The speed of sound in the medium remains close to that of an ideal gas, as shown in 
figs, l(e-f) and 2(e-f). As expected, in the massless gluon case, (? s = 1/3. The deviation 
of the dot-dashed curve from the ideal gas result is due to the finite current strange quark 
mass even though the light quarks and gluons are massless. In this case, the system moves 
further from the ideal gas behavior at later times as the temperature becomes comparable 
to m s . When the system is composed of massive gluons only, the deviation from an ideal 
gas is largest because the effective gluon mass is larger than the effective quark masses. 
The addition of the lighter massive quarks into the system tends to bring the sound speed 
closer to the ideal gas value. At the LHC, the higher temperature keeps the system closer 
to the ideal gas limit than at RHIC. 

The volume of the system, crucial to the determination of the charm yield, increases 
as shown in figs, l(g-h) and 2(g-h). Since the volume depends on the rapidity extent of 
the plasma, the relatively lower energy density and pressure of the gluon gas require a 
larger spatial extent to ensure that the energy of the system is equal to E tot . Because we 
have changed the initial conditions according to the composition of the minijet gas, the 
volume changes less than the parton gas volume. 

5 TOTAL CHARM YIELD 

In this section, we present our results on the cc production rates. The parton gas produces 
the largest cc yield because of the slower time evolution demonstrated in the previous 
section. Since the time evolution of the quark and gluon effective masses determines the 
relative enhancement of charm production by the massive excitations, we show the time 
dependence of the effective masses in a quark-gluon system in fig. 3. The finite current 
strange quark mass at zero temperature results in a slightly higher effective strange quark 
mass compared to the effective light quark mass. It also reduces the time dependence of 
the strange quark effective mass. The larger slope of the gluon effective mass as a function 
of temperature results in a faster decrease in gluon mass as a function of time. We have 
shown the effective masses for as long as T > T c - the minijet gas at RHIC is at T c when 
r ~ 1.5 fm. At later times, the finite masses become less effective for producing charm. 
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We remark that the running coupling constant, g, is also a function of time. The weaker 
coupling at later times also reduces the yield. 



The initial rate is approximately a factor of two larger when the quark degrees of 
freedom are included. At the beginning of the evolution, the rate is nearly independent 
of the initial parton mass. However, since the system cools more slowly with massive 
initial partons, the rate is larger at later times in the massive case. The number of cc 
pairs produced during the lifetime of the plasma is found by multiplying the rate, eq. 
(11), by the volume from eq. (30). In most cases, charm production only occurs during 
the early part of the evolution. At RHIC, production by the parton gas is essentially over 
after ~ 3 fm while production from the minijet gas is ended by ~ 0.5 fm. Note however, 
that cc pairs continue to be produced at much later times when the partons are massive, 
especially at the LHC. This is particularly true for the massive quark-gluon gas. 

The charm production rates and the number of produced cc pairs as a function of 
time at RHIC are given in fig. 4 for m c = 1.2 GeV and fig. 5 for m c = 1.5 GeV. The 
final charm pair yield in both cases is given in Table 1. The parton gas results are given 
in fig. 4(a-b) and 5(a-b). Here, although the production rate is larger when both quarks 
and gluons are included, the final number of cc pairs produced during the evolution of 
the system does not strongly depend on the composition of the plasma at RHIC because 
the larger rate is compensated by a correspondingly smaller volume. The yield from a 
massless quark-gluon gas is about 15% larger than that from a massless gluon gas. For 
massive initial quarks and gluons, the difference in the composition changes the yield by 
only 8%, as seen in Table 1. When m c = 1.2 GeV, the yield is increased 44% for a massive 
gluon gas relative to a massless gluon gas and 32% in the massive quark-gluon system. 
If m c = 1.5 GeV, the enhancement due to the massive partons is reduced to 36% for the 
gluons alone and to 29% for the quark-gluon gas. The yield at the lower mass is about 
3.5 times larger than when m c = 1.5 GeV. Some of the enhancement can be accounted 
for by the slower temperature evolution. After the system has evolved for 10 fm, the 
temperature of the massive gluon gas is approximately 10% larger than the massless gas. 
The difference is 8% when the quarks are included. Note particularly that the 5 thermal 
cc pairs produced with m c = 1.2 GeV is only a factor of two less than that expected from 
the initial production at this energy [13| . 



The yield is much smaller from a minijet gas due to the shorter equilibration time 
and the lower initial temperature. (In the minijet gas, To is reduced 24% for gluons alone 
and 53% for a quark-gluon system relative to the initial temperature of the parton gas.) 
In the minijet gas, the charm yield is reduced by a factor of 60-70 (for gluons) and 250-300 
(for quarks and gluons) compared to the parton gas. The thermal charm yield from the 
minijet gas is thus negligible compared to the initial charm rate [I3|. The influence of 



the massive partons is also reduced for the minijet gas. For m c = 1.2 GeV, the yield is 
only increased over the massless case by 23% in a gluon gas and 14% in a quark-gluon 
gas. The influence of the charm mass on the yield is also stronger for the minijet gas — the 
yield decreases by a factor of 4-5 with the larger charm mass. 
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The enhancement due to the effective parton mass is more substantial at the LHC, 
as seen in figs. 6 and 7 and in Table 2. The massive quarks and gluons have a more 
significant effect on the temperature evolution of the parton gas. The temperature is 27% 
higher after 10 fm with massive gluons and 20% higher for massive quarks and gluons. 
The enhancement of the yield is also larger: when m c = 1.2 GeV, the yield is 94% larger 
in the gluon gas and 70% higher in the quark-gluon gas. Increasing the charm mass 
only reduces the total yield by a factor of 2.5, thus the yield is less dependent on the 
charm mass at the higher temperature. Note that here the charm yield from the massive 
quark-gluon system, 250 pairs after 10 fm, is similar to the initial nucleon-nucleon rate 

0. 



The minijet gas is more effective at producing charm at LHC than at RHIC. Al- 
though the initial time remains short, To is 40% larger in the minijet gluon gas and 10% 
larger in the minijet quark-gluon gas than in the parton gas. Therefore the minijet charm 
yield is only a factor of two to six smaller than the parton gas yield and the minijet ther- 
mal charm yield also becomes a significant fraction of the initial production. Note also 
that the enhancement due to the massive partons is largest in the minijet gluon gas-a 
factor of 2.5 increase over that from a gas of massless gluons because the temperature of 
the massive gluon gas is 50% larger. 



In Ref. |13||, with an ideal, massless, quark-gluon plasma, one thermal charm pair 



was found at RHIC and 26 charm pairs at LHC with m c = 1.2 GeV. In that work, the 
initial temperature was nearly the same as the parton gas To used here but the initial 
time was closer to that of the minijet gas, r ~ (3T ) _1 . However, a larger, constant, 
g 2 kept the yield from being reduced. When initial conditions identical to those used 



here are chosen for the calculation of Ref. [13|, the results are quite similar. While the 



details of the expansion are somewhat different, we have checked that our massless charm 
production cross sections are in exact agreement with those in Ref. 



To study the dependence of the enhancement on the charm pair mass, we also 
calculate the thermal charm pair mass distributions for RHIC in fig. 8 and for the LHC in 
fig. 9. Generally, the thermal charm mass distributions are steeper than those charm pairs 
produced in the initial nucleon-nucleon interactions. (See ref. [[TB[] for the initial charm 
pair mass distributions.) The shapes of the mass distributions for each of the four cases 
we have studied with our two sets of initial conditions are quite similar, especially for the 
parton gas where T is the same in all cases. The shapes of the distributions from the 
minijet gas at RHIC are also similar although the initial yield is larger from the gluon gas 
because of its higher To. The enhancement is generally largest for low mass charm pairs, 
not far above threshold. At LHC, where the enhancement is greatest, the distributions 
with the massive excitations included approach the mass distributions for the massless 
case only at M ~ 8 - 10 GeV. 



12 



6 DISCUSSION 



We have investigated a new mechanism for enhancing thermal charm production by a 
quark-gluon plasma: massive excitations in the plasma state. We chose two different sets 
of initial conditions, a parton gas and a minijet gas, and calculated the thermal charm 
yield from each for both massless and massive quarks and gluons. In our calculation we 
assumed that the system stayed in thermal and chemical equilibrium during the expansion. 
Therefore, our charm yield is an upper bound on secondary charm production. 

The largest charm production was found when a parton gas was assumed because the 
characteristic thermalization time was r = 0.5 — 0.7 fm. Then the parton gas lifetime was 
long, more then 10 fm. At RHIC we obtained a 30-40% enhancement with massive gluons 
and quarks while at the LHC, a 50-100% enhancement may be expected. With m c = 1.2 
GeV, the charm mass used in recent calculations, we obtained 4.9 secondary charm quark 
pairs at RHIC and 245 charm pairs at LHC. Note that these numbers, upper limits on 
secondary charm production, are similar to the expected initial charm production If 
we take m c = 1.5 GeV, secondary charm pair production is reduced to 1.3 pairs at RHIC 
and 94 pairs at LHC. 

The lifetime of the minijet gas was significantly shorter than the parton gas due to 
the short thermalization time, r ~ 0.1 fm. This difference strongly reduced the charm 
yield from the minijet gas compared to the parton gas. The typical lifetime of the minijet 
gas was of order 2.5 fm although at RHIC the temperature of the quark-gluon gas dropped 
below T c after only 1.5 fm. The fast expansion reduces the influence of the massive quasi- 
particles at RHIC where the initial temperature of the minijet gas was also smaller than 
the initial temperature of the parton gas. The enhancement was typically 20% with 
massive quarks and gluons but the total yield, 0.016 pairs from the quark-gluon gas, was 
very small compared to the initial charm rate. Although the expansion was also fast 
at the LHC, the much higher initial temperature, T ~ 1 GeV, generated large charm 
production rates. With m c = 1.2 GeV we obtained 38 secondary charm quark pairs from 
a quark-gluon gas and 100 pairs from a gluon gas while with m c = 1.5 GeV, the yield 
was 16 and 47 pairs respectively. Note that the minijet gluon gas result is only a factor 
of two smaller than the corresponding parton gas yield at LHC. 

Our results show that, regardless of the initial conditions, the massive excitations of 
the quarks and gluons significantly enhance charm production by the plasma. Introducing 
thermal gluon and quark masses slowed the expansion of the system. The longer lifetime 
of the plasma as well as the reduced threshold for charm production with massive quarks 
and gluons leads to an enhancement of charm production over the massless case. Thus 
charm enhancement in heavy ion collisions could be an excellent probe of the presence 
of collective excitations in the deconfined plasma. If such enhanced charm production, 
beyond that predicted from the initial nucleon-nucleon collisions is observed, it could be 
expected that other exotic phenomena due to massive quark and gluon excitations may 
be found. 
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RHIC 




m c = 1.2 GeV 


m c = 1.5 GeV 




parton gas 


minijet gas 


parton gas 


minijet gas 


m = g 


3.2 


0.053 


0.93 


0.0128 


m/Oj 


4.6 


0.065 


1.27 


0.0154 


m = g + q 


3.7 


0.014 


1.07 


0.0027 


171 7^ g + q 


4.9 


0.016 


1.38 


0.0030 



Table 1: The total thermal cc pair yield from a parton gas and a minijet gas at RHIC. 
We consider both a massless and massive pure gluon gas and a quark-gluon system with 
massless and massive components. 
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LHC 




m c = 1.2 GeV 


m c = 1.5 GeV 




parton gas 


minijet gas 


parton gas 


minijet gas 


m = g 


102 


39 


43 


21 


m/Oj 


198 


101 


76 


47 


m = g + q 


145 


22 


60 


9.7 


171 7^ g + q 


245 


38 


94 


15.6 



Table 2: The total thermal cc pair yield from a parton gas and a minijet gas at LHC. 
We consider both a massless and massive pure gluon gas and a quark-gluon system with 
massless and massive components. 
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Figure Captions 



Fig. 1: The time evolution of the plasma produced in Au+Au collisions at RHIC is 
examined for a parton gas in (a), (c), (e) and (g) and a minijet gas in (b), (d), (f) 
and (h). The temperature is shown in (a) and (b), the energy density in (c) and 
(d), the square of the sound speed in (e) and (f) and the plasma volume in (g) and 
(h). The solid curve is for a massless gluon gas only while the dashed curve is for a 
gas of massive gluons. The dot-dashed curve is calculated assuming massless gluons 
and light quarks and m s = 150 MeV for the strange quark. The dotted curve is the 
result when the quarks and gluons have effective masses. 

Fig. 2: The same as fig. 1 for Pb+Pb collisions at LHC. 

Fig. 3: Time dependence of the effective quark and gluon masses. The parton gas results 
are given in (a) RHIC and (b) LHC. The minijet gas results are shown in (c) RHIC 
and (d) LHC. The RHIC results are given for Au+Au collisions, the LHC results 
for Pb+Pb collisions. The solid curve is the effective light quark mass, the dashed 
is the strange quark mass. The dot-dashed curve is the effective gluon mass. 

Fig. 4: The production rate (a), (c) and cc pair yield (b), (d) are given for Au+Au 
collisions at RHIC with m c = 1.2 GeV. The parton gas results are given in (a) and 
(b) while the minijet gas results are shown in (c) and (d). The solid curve is for 
a massless gluon gas only while the dashed curve is for a gas of massive gluons. 
The dot-dashed curve is calculated assuming massless gluons and light quarks and 
m s = 150 MeV for the strange quark. The dotted curve is the result when the 
quarks and gluons have effective masses. 

Fig. 5: The same as in fig. 4 for Au+Au collisions at RHIC with m c — 1.5 GeV. 

Fig. 6: The same as in fig. 4 for Pb+Pb collisions at LHC with m c = 1.2 GeV. 

Fig. 7: The same as in fig. 4 for Pb+Pb collisions at LHC with m c = 1.5 GeV. 

Fig. 8: The cc pair mass distribution from Au+Au collisions at RHIC with m c = 1.2 
GeV (a), (b) and m c — 1.5 GeV (c), (d). The parton gas results are given in (a) 
and (c) while the minijet gas results are shown in (b) and (d). The solid curve is 
for a massless gluon gas only while the dashed curve is for a gas of massive gluons. 
The dot-dashed curve is calculated assuming massless gluons and light quarks and 
m s = 150 MeV for the strange quark. The dotted curve is the result when the 
quarks and gluons have effective masses. 

Fig. 9: The same as in fig. 8 for Pb+Pb collisions at LHC. 
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